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5^ ' The hadronic fc^ -spectrum inside a high energy jet is determined including corrections of relative 

^ I magnitude O (y^ol) with respect to the Modified Leading Logarithmic Approximation (MLLA), in 

the limiting spectrum approximation (assuming an infrared cut-off Qo = ^qc^ ) s-nd beyond (Qo 7^ 
^qcd)' '^^^ results in the limiting spectrum approximation are found to be, after normalization, in 
impressive agreement with preliminary measurements by the CDF collaboration, unlike what occurs 
at MLLA, pointing out small overall non-perturbative contributions. Within the same framework, 
2-particle correlations inside a jet are also predicted at NMLLA and compared to previous MLLA 
f~i ■ calculations. 
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INTRODUCTION 



(N 
> 

The production of jets - a coUimated bunch of hadrons - in e+e , e p and hadronic colhsions is an ideal playground 
to investigate the parton evolution process in perturbative QCD (pQCD). One of the great successes of pQCD is the 
existence of the hump-backed shape of inclusive spectra, predicted in [1] within the Modified Leading Logarithmic 
Approximation (MLLA), and later discovered experimentally (for review, see e.g. [2]). Refining the comparison of 
pQCD calculations with jet data taken at LEP, Tevatron and LHC will ultimately allow for a crucial test of the Local 

I Parton Hadron Duality (LPHD) hypothesis [3] and for a better understanding of color neutralization processes. 

^~») Progress towards this goal has been achieved recently. On the theory side, the inclusive fc^^ -distribution of particles 
• • . inside a jet has been computed at MLLA accuracy [4], as well as correlations between two particles in a jet [5]. Analytic 
^ calculations have first been done in the limiting spectrum approximation, i.e. assuming an infrared cutoff Qo equal 
to A^^^ (A = ln(5o/A^^^ = 0). Subsequently, analytic approximations for correlations were obtained beyond the 
5—1 ] limiting spectrum using the steepest descent method [6] . Experimentally, the CDF collaboration at Tevatron reported 
on fc^ -distributions of unidentified hadrons in jets produced in pp collisions at ^/s = 1.96 TeV [7]. 

MLLA corrections, of relative magnitude O {^/c^T) with respect to the leading double logarithmic approxima- 
tion (DLA), were shown to be quite substantial for single-inclusive distributions and 2-particle correlations [4, 5]. 
Therefore, it appears legitimate to wonder whether corrections of order 0{as), that is next-to-next-to- leading or 
next-to-MLLA (NMLLA), are negligible or not. 

The starting point of this analysis is the MLLA evolution equation for the generating functional of QCD jets [8]. 
Together with the initial condition at threshold, it determines jet properties at all energies. At high energies one can 
represent the solution as an expansion in y^aj. Then, the leading (DLA) and next-to-leading (MLLA) approximations 
are complete. The next terms (NMLLA) are not complete but they include an important contribution which takes 
into account energy conservation and an improved behavior near threshold. An example of a solution for the single 
inclusive spectrum from the MLLA equation is the so-called "limiting spectrum" (for a review, see [8]) which represents 
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a perturbative computation of the spectrum at A = with complete leading and next-to-leading asymptotics. Some 
results for such NMLLA terms have been studied previously for global observables and have been found to better 
account for recoil effects. They were shown to drastically affect multiplicities and particle correlations in jets: this 
is in particular the case in [9], which deals with multiplicity correlators of order 2, and in [10], where multiplicity 
correlators involving a higher number of partons are studied; in particular, the higher this number, the larger turn 
out to be NMLLA corrections. 

The present study makes use of this evolution equation to estimate NMLLA contributions to our differential 
observables. It presents the complete calculations of the single inclusive fc^-distribution leading to the main results 
published in [11], and extends them to 2-particle correlations inside a high energy jet. 

The paper is organized as follows. First, Section II presents a system of evolution equations including O (ag) 
corrections, which allows for the computation of the inclusive spectrum, G, beyond MLLA accTiracy. Section III 
is devoted to the NMLLA evaluation of the color currents of quark and gluon jets and, from them, to the inclusive 
A;_L-distribution in the limiting spectrum approximation. These predictions are also compared to preliminary measure- 
ments performed recently by the CDF collaboration. Going beyond the limiting spectrum is the subject of Section IV, 
in which inclusive fcj^ -distributions are computed at an arbitrary A. The 2-particle correlations including NMLLA 
corrections are determined in Section V. Finally, the present approach and the results obtained in this paper are 
discussed in detail and summarized in Section VI. 



II. EVOLUTION EQUATIONS 



A. Logic and energy conservation 



As a consequence of the probabilistic shower picture, the notion of Generating Functional (GF) was proved suitable 

to understand and include higher order corrections to DLA asymptotics (sec [8] and references therein). 

The single inclusive spectrum and the n-particle momentum correlations can be derived from the MLLA Master 
Equation for the GF Z = Z{u) [8] after successively differentiating with respect to a certain probing function u = u{k); 
k denotes the quadri-momentum of one parton inside the shower and the solution of the equations are written as 
a perturbative expansion in as- At high energies this expansion can be resummed and the leading contribution be 
represented as an exponential of the anomalous dimension 7(05); since further details to this logic can be found in 
[5, 8], we only give the symbolic structure of the equation for the GF and its solution as 

^^My)Z Z^exp|y%(a,(y')W} (1) 

where j{as) can be written as an expansion in powers of ^/a^ 

7(a«) = ^^+as+ + + . . . (2) 

The equation in (1) applies to each vertex of the cascade and its solution represents the fact that successive and 
independent partonic splittings inside the shower, such as the one displayed in Fig. 1, exponentiate with respect to 
the evolution-time parameter dy = dQ/Q; 9 <C 1 is the angle between outgoing couples of partons. The choice of 
y follows from Angular Ordering (AO) in intrajet cascades; it is indeed the suited variable for describing time-like 
evolution in jets. Thus, Eq. (1) incorporates the Markov chains of sequential angular ordered partonic decays which 
are singular in © and 7(0:5) determines the rate of inclusive quantities growth with energy. 

While DLA treats the emission of both particles as independent by keeping track of the first term ^ ^/ci^ in (2) 
without constraint, the exact solution of the MLLA evolution equation (partially) fulfills the energy conservation in 

each individual splitting process (^; -|- (1 — 2;) = 1) by incorporating higher order (a"^^, n > 1) terms to the anomalous 
dimension. Symbolically; the first two analytical steps towards a better account of these corrections in the MLLA, 
NMLLA evolution, which we further discuss in II C, can be represented in the form 



A7 ~ y {as + asi~'^ \nz)dz ~ -h a^/^. 



where i = ln(l/a;) ~ as with a; ^ 1 (fraction of the jet energy taken away by one hadron), 2; ~ 1 for hard partons 
splittings such as g ^ qq. . . (this is in fact the region where the two partons are strongly correlated). 

Energy conservation is particularly important for energetic particles as the remaining phase space is then very lim- 
ited. On the other hand, a soft particle can be emitted with little impact on energy conservation. Some consequences 
of this behavior have also been noted in [12]: 
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(i) the soft particles follow the features expected from DLA; 
(ii) there is no energy dependence of the soft spectrum; 

(iii) the ratio of soft particles r = Ng/Nq in gluon and quark jets is consistent with the DLA prediction Nc/Cp = 9/4 
(see the measurement by DELPHI [13]). 

This is quite different from the ratio of global multiplicities which acquires large corrections beyond DLA (see, 
for example Fig. 18 in the second reference given in [14]). For this quantity the HERWIG parton shower model 
corresponding to MLLA and exact energy conservation (same Fig. 18) and the full summation of the perturbative 
series of MLLA evolution equation (see also [15]) come close to the data at r = Ng/Nq « 1.5 at LEP energies. As an 
intermediate example, we can mention the successful description of the semi-soft particle ln(l/x) distribution ("hump- 
backed plateau") where the first correction (MLLA), despite the large value of the expansion parameter ^/a^ « 0.35, 
already gives a good description of the data at the peak {Q = 91.2 GeV) of the e+e~ annihilation into a qq pair 
[16]. 



B. MLLA evolution 



We study the formation of hadrons inside a jet produced in high-energy scattering processes, such as e+e~ anni- 
hilation or pp and pp collisions. A jet of total opening angle 8o is initiated by a parton A (either a quark, Q, or a 
gluon, G) with energy E; A then splits into partons B and C, with energy fractions z and (1 — z) respectively, forming 
a relative angle 6 (see Fig. 1). At the end of the cascading process, the parton B fragments into a hadron h with 
energy xE, with the fragmentation function 

B{z) = -D^ (-,zEeo, Qo) , iB = Q, G) (3) 

Z \Z I 

which describes the distribution of the hadron h inside the sub-jet B with an energy-fraction xjz. As a consequence 



h 




Figure 1: Parton A with energy E splits into parton B (respectively, C) with energy zE (respectively, (1 — «)-£) which fragments 
into a hadron h with energy xE. 

of AO in parton cascades, the functions Q(z) and G(z) satisfy the system of two-coupled integro-differential evolution 
equations [5]: 



(4) 



$9( 
9^ 



Qy — 


dQ 
dy 


[ dz 
Jo 


as 

IT 


Gy = 


dG 
dy 


[ dz 
Jo 


as 

TT 



;(z)(1-z)(g(z) + G(1-z)-g) 
+ nf ^liz) {2Q{z)-G) 

with as, the running coupling constant of QCD, given by 

as = as{i,y) = 



27r 



4N,f3o{e + y + X)' 



(5) 



(6) 
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and where we define 

^ = ln(l/.), y = ln^=ln^, A = ln-^, (7) 

following the notations of Ref. [8] ; the MLLA equations above follow from the GF logic commented in the introductory 
paragraph. The scale Qo appearing in (7) is the coUinear cut-off parameter, A is the non-perturbative scale of 
QCD which we set to 250 MeV in this work [32] , and 



is the first term in the perturbative expansion of the /3- function (A^^ is the number of colors, Tr = n//2 where n/ = 3 
is the number of light quark flavors). We only consider in this work the 1-loop expression for the running coupling 
constant, assuming that the role of the conservation of energy is much more important than the effects of 2-loop 
corrections to a^, as seen for instance in the case of multiplicity distributions [14]; we shall discuss this further in 
Section VI A. The coupling constant as is also linked to the DLA anomalous dimension 70 of twist-2 operators by 

,lil.y) = 2N.^ = Ye = i + y = In^. (9) 

In Eqs. (4) and (5), ^^{z) represent the one-loop DGLAP splitting functions [8] and wc note: 

Q = Q(l) = xDi^ix, EQo, Qo), G = G(l) = a;Z)^(a;, EQo,Qo)- 
In the small a; -C -2 limit which we consider here, the fragmentation functions behave as 



B{z) « pMln-,ln— ^ = p>^ {\nz + i,y) , (10) 
\ X Qq J 

where pg is a slowly varying function of the two logarithmic variables \n.{z/x) and y [1] that describes the hump-backed 
plateau. 



C. Taylor expansion 



The resummation scheme at MLLA is discussed in [5], in which G{z) and G{1 — z) were replaced by G(l) in the 
non-singular part of the integrands in Eqs. (4) and (5). In the present work, we calculate next-to-MLLA (NMLLA) 
corrections from the Taylor expansion of pg in the variables Inz and ln(l — z) in the domain: 



corresponding to hard parton splittings. To first order, 

dpiln z) 



p(ln z) = p(ln z = 0) + 



dlnz 



In 2=0 



\ lnz\ ~ |ln(l-z) 
lnz + 0{\n^ z) , 



p(ln(l-z)) = p(ln(l-z) = 0) 
or, equivalently, for the function B{z): 



9p(ln(l - z)) 



d\n{l-z) 



ln(l-z)=0 



\n{l-z) + 0{\n^{l-z)) 



lnz|<^ 



B{z) ' " B{l)+Be{l)\nz + 0{ln^ z) 

|ln(l-z)|<^ 



B{1 - z) 



B{1) + Be{l) ln(l -z) + {ln\l - z)) . 



(11) 

(12) 

(13) 
(14) 



The derivative with respect to Inz or ln(l — z) has been replaced by the one with respect to i because of (10) and 
the property that, at low x, B is a. function of {Inz + 1) or (ln(l — z) -\- 1). Since £ = O (l/y^) (see [8]) the above 
expansion can be written symbolically 



B (z) - B (1 - z) ~ ci C2(VSr) + 0{as 



Cl, C2 



0{l). 



The terms proportional to Bi thus provide NMLLA corrections to the solutions of the MLLA evolution equations (4) 
and (5). 
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D. Evolution equations including NMLLA corrections 

1. Quark jet 

In order to determine NMLLA corrections to the evolution equation (4), the 1-loop splitting functions (see [8]) are 
written 

^l{z) = Cf{^- + , (1 - z)^l{z) = 27V, Q + 4>liz?j , 

where (jy^iz) = {z — 2) and (j)^{z) = {z — 1) {2 — z{l — z)) are regular functions of z. The term proportional to G{z) 
in the integrand of (4) becomes 



d^ ^ mz) G{z) = 2Cf f—'^ G{z) + Cf f dz ^ cj^^{z) G{z), 

Jo ^ Jo Z TV Jo ^ 

the second part of which is expanded according to (13). Replacing as/n = Jq/2Nc (see 9), one gets 



I dz^niz)Giz) . f [(/ '-l,^,Giz))-l,^G^l,^,G,^... 



(15) 



(16) 



where Gi = G£{1) and Qe = Qe{l)- The first integral in the r.h.s of (16) provides the DLA (leading) term as 
2; — > 0, while the second and third terms correspond to higher powers of ^/a^, that is MLLA and NMLLA corrections 
respectively. The z-dependence of ag in (16) has only been taken into account in the singular (DLA) part dominated by 
small z. On the contrary, for the non-singular parts corresponding to branching processes in which z ^ 1 — z — O (1), 
as has been taken out of the z integral [33] as done in [5]. The dependence on the other variables, kj_, 0, is of course 
unchanged. 

Likewise, the term proportional to Q{1 — z) — Q in (4) can be expanded according to (13), leading to 



/' dz ^ mz) (q{1 -z)-Q) ^ f dz ^Q, ^{z) ln(l 
Jo ^ ^ ' Jo ^ 



Cf 

Nr. 



Gi. 



(17) 



In the second line of (17), we have used the approximated formula Qi w Cf/Nc G^-|-0(7q) that holds at DLA because 
subleading terms would give C(7o) corrections which are beyond NMLLA (see also appendix A). Finally, plugging 
(16) and (17) into (4), we obtain 



7 CV/5_7r2 

8 iVc I 8 ~ Y 



(18) 



where the term proportional to ^^Gg = 0(70) constitutes the new NMLLA correction. It is quite sizable and should 
be taken into account in the coming calculations. 



2. Gluon jet 



Along similar steps, we now evaluate NMLLA corrections to Eq. (5). The first term in the integral can be cast in 
the form 

/ dz^ $9(z)(1-z)(g(z) + G(1-z)-G) 

°/^,j„„)_ii,jo.(g4).jG„ («, 

«, [ d. "-^ (2«(,) - g) . -^,1 (2Q - G) - I'^llQ,. (20) 



and the second into 
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Summing (19) and (20), replacing like before Q by its DLA formula Q s 
details), the evolution equation for particle spectra inside a gluon jet reads 



Cp/Nc G (see appendix A for further 



Gy — 



z 



"11 


nfTn 




.12 + 


3iVe 


/67 




13 


V 36 


6 


18 



Nr. 



llG 



Nr. 



(21) 



The first term in parenthesis in (18) and (21) is, as stressed before, the main (double logarithmic) contribution. 
According to the Low-Barnett-Kroll theorem [17], the dz/z term, which is of classical origin, is universal, that is, 
independent of the process and of the partonic quantum numbers. The other two (single logarithmic) contributions, 
which arise from hard parton splitting, are quantum corrections. It should also be noticed that, despite the large size 
of NMLLA corrections coming from g ^ gg and g ^ qq splittings, a large cancellation occurs in their sum (21). The 
coefficients of the terms proportional to Gi in (18) and in (21) are in agreement with [18]. 



3. NMLLA system of evolution equations 

Once written in terms of £' = ln{z/x) and y' = \n{xEQ/Qo), the system of two-coupled evolution equations (18) 
and (21) finally reads, 



Q{£, y) = 6{i) + ^ t di' r dy'-f^ii' + y') fl - ai<5(f - £) + h^Sil' - y')] G(£', y'), 



G{l,y) 



rv 

5{l) + di' dy'- 
Jo Jo 



7o'(^' + y') l-a^S{e'-i)+a2S{e'-e)M^',y') G{£',y'), 



with ipe = Gi/G and the MLLA and NMLLA coefficients [34] given by: 



ai 

0.2 

a2 



3 

7 Cf (h 7r2\ 

- + — « 0.42, 

8 TVc V8 6 y 

67 TT^ 13n/TiiCi. n/=3 



36 



18 Nr. Nr 



0.06. 



(22) 

(23) 

(24a) 
(24b) 

(24c) 

(24d) 



As can be seen, the NMLLA coefficient 02 is very small This may explain a posteriori why the MLLA "hump-backed 
plateau" agrees very well with experimental data [1, 19]. Therefore, the NMLLA solution of (23) can be approximated 
by the MLLA solution of G {i.e. taking 02 = 0), which will be used in the following to compute the inclusive kj_- 
distribution as well as two-particle correlations inside a jet [35]. The MLLA gluon inclusive spectrum is given by 
[8]: 



Gii,y)=2 



TiB) p dr 
/3i 



B) n dr 
'0 Jo 



^B{T,y,i), 



(25) 



1 y 

where the integration is performed with respect to r defined by a = - In - + ir and with 



:FB{T,y,e) = 



cosh a 



sinha 



£ + y a 
Po sinh a 



B/2 



/s (2 V^(t, ?/,£)), 



7 



Z{T,y,£) 



i + y a 
Po sinh a 



( y- (■ \ 
cosh a sinh a 1 

V V + ^ J 



B = ai/Po and Ib is the modified Bessel function of the first kind. To get a quantitative idea on the difference 
between MLLA and NMLLA ghion inclusive spectrum, the reader is reported to Appendix B where a simphfied 
NMLLA equation (23) with a frozen couphng constant is solved. The magnitude of 02, however, indicates that the 
NMLLA corrections to the inclusive quark jet spectrum may not be negligible and should be taken into account. After 
solving (23), the solution of (22) reads 



Q{i, y) = ^ y) + («i - ai)Gi{l, y) + (ai (ai - ai) + - aa) Gu{l, y)] + O(7o) 



(26) 



It differs from the MLLA expression given in [4] by the term proportional to Gu, which can be deduced from the 
subtraction of (Cf/AT^) x (23) to Eq. (22). 



III. SINGLE-INCLUSIVE A; -DISTRIBUTION IN THE LIMITING SPECTRUM 

While MLLA calculations show that, asymptotically, the shape of the inclusive spectrum becomes independent of 
A [8, 20], sotting the infrared cutoff Qo of cascading processes as low as the intrinsic QCD scale A^^^ is a daring 
hypothesis, since it is tantamount to assuming that a perturbative treatment can be trusted in regions of large running 
as- However, it turns out that, experimentally, this shape is very well described by A = 0. We shall show below that 
this remarkable property is also true for the single-inclusive fc_L-distribution. This will be further confirmed in section 
IV in which non-vanishing values of A are considered. 



A. Double-differential distribution 



The double differential distribution d^A^/(dxdln6') for the production of a single hadron h at angle 8 in a high 
energy jet of total energy E and opening angle > G, carrying the energy fraction x, is obtained by integrating the 
inclusive double differential 2-particle cross section (see [4]) [36]. Then, the single- inclusive fc_L-distribution of hadrons 
inside a jet is obtained by integrating d^A''/(dxdln^) over all energy-fractions x: 



din k± 



As in [4], a lower bound of integration, iraim is introduced since the present calculation is only valid in the small-a; 
region, and therefore cannot be trusted when £ = ln(l/a;) becomes "too" small. We shall discuss this in more detail 
in Section III C and Appendix G. 

According to [4] , jJ^j^q can be expressed as 



d^N d 



da;dlne dlnG 



Fl:^{x,e,E,eo), (28) 



where F/^, which represents the inclusive production of h in the sub-jet of opening angle 6 inside the jet Aq of opening 
angle 60, is given by a convolution product of two fragmentation functions [4]: 



Fj^^ {x, e, E, Qo)=J2 [ du Dl {u, EQq, uEQ) D'I (-,uEe, Q, 



(29) 



The convolution expresses the correlation between the energy flux of the jet and one particle within it. Eq. (29) is 
schematically depicted in Fig. 2: u is the energy-fraction of the intermediate parton A, D'^ describes the probability to 
emit A with energy uE off the parton Aq (which initiates the jet) taking into account the evolution of the jet between 
Qq and 0, and describes the probability to produce the hadron h off A with energy fraction x/u and transverse 
momentum k±_ « uEQ > Qo; k±_ is defined with respect to the jet axis which is, in this context, identified with the 
direction of the energy flux. 



Figure 2: Inclusive production of hadron h at an angle Q inside a high energy jet of total opening angle Qq. 



As discussed in [4], the convolution (29) is dominated by u = O {!). Therefore, D^^ {u, EQo,uEQ) is given by 
DGLAP evolution equations. On the contrary, since x <^ u = O {1) in the small-a; limit where MLLA evolution 
equations are valid, behaves as (see (10)) 



Dl (^,uEe, Qo) %1 (in ^, Inu + Ye) 



(30) 



x<« U 
' J X^ " \ X 

Since le + Inu = ^ + Inu + y, the hump-backed plateau p\ depends on the two variables ^ + Inu and y, and we 



conveniently define D as: 



u \u / 



(31) 



The Taylor expansion of p'l to the second order in \nu for m~ 1<^ |lnM| <C 1, that is, one step further than in [4], 
leads to 



xFl^{x,Q,E, Go) « xFl^{x,e,E,eo) 



IE 



du u{lnu)Dl{u, EGo, uEQ) 



d^2 



where 



xF^{x,e,E,Qo)^Yl [ duu(l + ilnu)i;,4e,y))Dl{u,Eeo,uEe) 



is the MLLA distribution calculated in [4]. In (33) we have introduced first logarithmic derivatives of D'l 

1 db,{t,y) 



Dl{l,y) dl 



(32) 



(33) 



(34) 



Thus, as in [4], in the soft limit the correlation disappears and the convolution (29) is reduced to the factorized 
expression in (32). 

The second term in the r.h.s. of (32) is the new NMLLA correction calculated in this paper. Since x/u is small, the 
inclusive spectrum D'l^t, y) occurring in (32) should be taken as the next-to-MLLA solution of the evolution equations 
(22) and (23). However, as already mentioned and shown in Appendix B, the MLLA inclusive spectrum for a gluon 
jet can be used as a good approximation for (23) (with oi 7^ 0, n2 = 0) such that, in (33), it is enough to use this level 
of approximation. So, we shall therefore use Eqs. (25) and (26) in the following. 

The NMLLA correction in (32) globally decreases Ixi^jJ'l in the perturbative region {y > 1.5). Indeed, while the 
MLLA part proportional to In?! in (33) is negative [4], it is instead, there, positive because of the positivity of u and 



the MLLA corrections when y is large enough, 



In^u and "j^" ~ ^p- (see Fig. 18 in Appendix C). The NMLLA contribution therefore tempers somehow the size of 



B. Color currents 



The function Fl^^ is related to the inclusive gluon distribution via the color currents defined as [4, 8] 

xF'^=^^G{l,y). (35) 
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The color current can be seen as the average color charge carried by the parton A due to the DGLAP evolution from 
Ao to A. Introducing the first and second logarithmic derivatives of D^, 



(V'M + V'A,«)(^,y) 



1 



which are MLLA and NMLLA corrections, respectively, Eq. (32) can now be written 



(36) 



where 



{u 



In* u)i^ = I du{u In* u) {u, EQq, uEQ) . 
Jo 



(37) 



Unlike in [4] at MLLA, using the approximation u = 0{1) to replace in (37) uEQ by EQ requires here some care, 

since the resulting scaling violation of the DGLAP fragmentation functions also provides 0{as) corrections to (u). 
Explicit calculations (see Appendix D) show that they never exceed 5% of the leading term. Accordingly, we neglect 
them in the following and replace (37) by 



L 



{u In' u}l ~ / du{u In* u) Dt {u, EQq, EQ) . 



(38) 



The total average color current {C)^^ of partons caught by the calorimeter decomposes accordingly into three terms 
which can be written: 



\MLLA-LO 



I r//0\lNJ 



NMLLA-MLLA 



(39) 



The leading order (LO) O (1) and MLLA O {^/a^) contributions to the color currents have been determined in [4]. 
The new NMLLA O (as) correction evaluated in this paper reads 



f //0\JNJ 



NMLLA-MLLA 



= iVe {uln^ u)l ii^le + i^gM) + Cf (win' u)1 {i^l^ + Vg,«), 



(40) 



assuming Q = Cp/Nc G. We checked that using instead the NMLLA exact formula (26) for the quark inclusive 
spectrum Q actually leads to negligible corrections to the color currents (see Appendix E). Eq. (40) can be obtained 
from the Mellin-transformed DGLAP fragmentation functions 



Jo 



through the formula 



J=2 



/ duuln^ uDI^{u,^). 

Jo 



(41) 



for quark and gluon jets are 



Given the rather lengthy expressions, the complete analytic results for ^(^^nmlla-mlla 
given in Appendix F. 

For illustrative purposes, the color currents are plotted in Fig. 3 in the limiting spectrum approximation (A = 0). 
The LO (solid line), MLLA (dash-dotted) and NMLLA (dashed) currents are computed for a quark (left) and for a 
gluon jet (right) with energy Yq^ = 6.4 - corresponding to Tevatron energies - and at fixed 1 = 2. As can be seen 
in Fig. 3, NMLLA O {ug) corrections to the MLLA color currents are clearly not negligible, yet of course somewhat 
smaller than the MLLA O (\/al) corrections to the LO result. In the pcrturbativc region {y > 1.5), these corrections 
are positive and consequently decrease the difference with the LO estimate. On the contrary, at small y < 1.5, the 
corrections are rather large and negative coming from the negative sign of Gu{(-,y) (see Fig. 18 in Appendix C). 
However, it should be kept in mind that as y goes to 0, kj_ gets closer to A^^^ (remind that Qq = A^^^ in the 

limiting spectrum approximation) and, thus, the present perturbative predictions may not be reliable in this domain. 

Note also that both the MLLA and NMLLA corrections vanish at y = (since Ge = Gu = 0) and when = Oq. 
Another interesting property to mention is the decrease of MLLA and NMLLA corrections as ^ increases, that is, 
when partons get softer and recoil effects more negligible. 
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From the color currents, the NMLLA double-differential 1-particle distribution at small x (see Eq. (28)) 

<1^N\ 1 ,„, d „,„ . 1 „,„ , d 



can be determined for any value of A. The NMLLA behavior of d?N /dlAy is therefore easily deduced from (C)ao and 
its y-dependence, d(C)Ao/dy. 



C. fcj^ -distributions 

The fcj^ -distributions of hadrons are computed from the numerical integration of the double-differential cross section, 
Eq. (42). On Fig. 4 are shown the MLLA (dashed lines) and NMLLA (solid lines) dN/dy distributions for a quark 
(left) and a gluon jet (right) with Ye^ = 4.3 and Yq^ = 6.4. The size of NMLLA corrections proves quite substantial 
over the whole y-range. We find in particular that at large y (or fc^^), the distributions at NMLLA are lower than 
at MLLA (and larger at small y). This softening of the spectra can be understood physically by the role of energy 
conservation in jets. With respect to DLA, MLLA and NMLLA take better into account the recoil of the emitting 
parton at each step of the cascading process. The fraction of energy carried away by the emitted soft partons gets 
reduced, which finally damps the final emission of hadrons at large fc_L [37]. As already stressed in Section III A, the 
value of the lower limit of integration £„iin below which the present small- a; calculation may not be trusted cannot be 
directly predicted. In [4], the appearance of positivity problems in the double-differential distribution at small t led us 
to consider a minimal value £min such that d^iV/d£dlnfc^ is kept positive for all i > £min, leading to [38] £niin — 2.5. 
For consistency, the same criterion is used in the present paper. We find that smaller values of € actually fulfill the 
positivity requirement, roughly iinm — 2 and im\n — 1 for quark and gluon jets at Tevatron energies. 

It is interesting to note that the range over which NMLLA calculations appear sensible extends to smaller £, 
therefore to larger x, than at MLLA; this also corresponds to larger y at fixed Y . One could therefore expect the 
present NMLLA predictions to agree with experimental results in a larger domain of k_\_. This is discussed in the 
coming Section. 



D. Comparison with CDF preliminary data 

The CDF collaboration at Tevatron recently reported on preliminary data of hadronic single-inclusive k_^- 
distributions inside jets produced in pp collisions at -/s = 1.96 TeV [7]. The measurements cover a wide domain 
of jet energies, with hardness Q — EQq ranging from Q — 19 GeV to Q = 155 GeV. The CDF results, including sys- 
tematic errors, are plotted in Fig. 5 together with the MLLA predictions of [4] (dashed lines) and the present NMLLA 
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Figure 4: MLLA (green) and NMLLA (blue) inclusive y-distributions for a quark (left) and a gluon jet (right) with Ybq = 4.3 
(Q = 19 GeV) and Fgo = 6.4 (Q = 155 GeV). 



calculations (solid lines). Data and theory are normalized to the same bin, Infc^ = —0.1, because of presently too 
large normalization errors in the CDF preliminary data. The experimental measurements reflect a mixing of quark 
and gluon jets: 

f dN \ f dN \ , J dN \ 

characterized by one Q-dependent mixing-parameter oj, estimated from PYTHIA [39], used in the theoretical calcu- 
lation. The agreement between the CDF results and the NMLLA distributions over the whole fc^ -range is excellent. 
The NMLLA calculation is in particular able to capture the shape of CDF spectra at all Q. Conversely, predictions 
at MLLA prove only reliable at not too large . 

The domain of validity of the predictions has been enlarged to larger (and thus to larger x since Y is fixed) 
computing from MLLA to NMLLA accuracy [40] . It is however to be mentioned that, due to the normalization at 
the first bin, this extension of the domain of prediction only concerns, strictly speaking, the shape of the distribution. 
Equally importantly, the agreement between NMLLA calculations and experimental results brings further support 
to the Local Parton Hadron Duality (LPHD) picture [3]. We indeed find it remarkable to observe that the entire 
fcj^ -domain probed experimentally can be very well described by strict perturbation theory, leaving out only limited 
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12 O 1 2 3 

In (k^ / 1 GeV) In (k^ / 1 GeV) 

Figure 5: CDF preliminary results on hadronic single-inclusive fe^ -distributions, compared with MLLA (dashed lines) and 
NMLLA (solid lines) calculations at the limiting spectrum; the boxes are the systematic errors (their lower limits are cut at 
large k± for the sake of clarity) 



non-perturbative dynamics in the production of hadrons inside a jet, at least for inclusive enough observable hke 
single-particle -distributions. 

E. Theoretical uncertainties 

The spectacular agreement between our NMLLA calculations and preliminary data should not hide the theoretical 
uncertainties that affect the former. 

First, we did not take into account all NMLLA corrections. While scaling violations have already been dealt with 
in subsection IllB and Appendix D, other NMLLA corrections arise from varying A^^^ in the expression of as- In 
Figs. 6 is plotted the inclusive fc_L-distribution {Q = 119 GeV) at values of A^^^ ranging from 150 to 500 MeV (left), 
as well as the ratio to its value at the default A = 250 MeV (right). All curves have been normalized to the bin 

QOD \ t> / 

ln(A:j^/lGeV) = —0.1. In the largest bin ln(fcj_/lGeV) = 3, varying A^^^ varies from 150 to 400 MeV does not yield 
a relative variation larger than 20%. The corresponding curves still fall within the experimental systematic errors. 
Note that the fact that variations seem only important at large k± only comes from the normalization procedure in 
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the bin ln(fc^/lGcV) — —0.1. A more delicate matter concerns the dominance of the type of NMLLA corrections 
that we have taken into account. Some remarks will concerning this point are postponed to the general discussion 
in section VI. The second point concerns the jet axis, which is defined here as the direction of the energy flow. It 




Figure 6: The dependence on A , absolute (left) and relative (right). 



is implicitly determined by a summation over all secondary hadrons in energy-energy correlations. At the opposite, 
the jet axis is experimentally determined exclusively from all particles inside the jet. Whether these two definitions 
match within NMLLA accuracy, O (as), is a matter which deserves further investigation. This goes however beyond 
the scope of the present work. 

Last, cutting the integral (27) at small I may look somewhat arbitrary. However, at the end of Appendix G, we 
provide in Figs. 21 curves which show the variation of the inclusive fcx-distribution at MLLA and NMLLA when 
^min is changed. Varying it from 1 to to 1.75 does not modify the NMLLA spectrum at large fcj^ by more than 20%. 
Variations are more dramatic at MLLA. 



IV. SINGLE-INCLUSIVE fc^ -DISTRIBUTIONS BEYOND THE LIMITING SPECTRUM 

A. Inclusive spectrum 

So far, the calculations have been performed in the limiting spectrum approximation, Qo = ■'^qcd or A = 0. This 
assumption, which cuts off hadronic yield below Qo should be valid as long as the mass of the produced hadrons is 
not too large as compared to A^^^ . This is the case when dealing mostly with pious. We perform in this Section the 
exact calculation of single- inclusive spectra as well as /c^ -distributions beyond this approximation, A 7^ 0, that is for 
hadrons with mass mn — Qo 7^ A ^ [20]. 

The inclusive gluon spectrum was given in [5] a compact Mellin representation: 

. ,,,,, / - ... (^)--(^)-\-. 
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from which an analytic approximated expression was found using the steepest descent method [6]. However, G{i,y) 
is here determined exactly from an equivalent representation in terms of a single Mellin transform (which reduces 
to (25) as A ^ 0) [20] 



G{e,y) = 



e + y + x 



e+zoo 



duj 



f3o B 2iTi 
X $(-A + B + 1, S + 2, -Lj{l + y + \)) IC{uj, A) 

which is better suited for numerical studies. The function /C appearing in Eq. (44) reads 



IC{lu,X) = 



r{B) 



{lu A)^ ^'{A,B + 1,lu A), 



(44) 



(45) 



where A = l/(/3o B = ui/Pq, and <i> and ^' are the confluent hypergeometric function of the first and second kind, 
respectively. The single-inclusive spectrum at MLLA is plotted in Fig. 7 for various values of A, A = 0, 0.2, 0.5, 1, for 
a gluon jet with Yq = 6.4. Increasing A reduces the emission in the infrared region and therefore favors hard particle 
production at ^ <C Y/2 (asymptotic position of the peak of the hump-backed plateau). Still, it is worth remarking 
that the global shape of G at finite A remains similar to that obtained in the limiting spectrum approximation. Note 
also that there is a discrete part at finite A, proportional to S{£), corresponding to the finite probability for no parton 
emission when Qq ^ A , the parton multiplicity becoming infrared finite at A 7^ (see the second reference in [20] ) . 



X 




Figure 7: Inclusive spectrum for a gluon jet (Yeo = 6.4) for different values of A. 



B. Color currents 



The color currents, Eq. (39), can now be determined beyond the limiting spectrum from the inclusive spec- 
trum calculated in the previous section. In Fig. 8 are displayed the MLLA corrections to the LO color current, 
^^^^MLLA-LO/^^^LO (jgf^)^ g^^^ NMLLA corrections to the MLLA color currents, j(c)NMLLA-mlla^^^^mlla (^^[g^t), 
for different values A = 0, 0.5, 1. Fig. 8 clearly indicates that the larger the values of A, the smaller the MLLA (and 
NMLLA) corrections. In particular, MLLA (NMLLA) corrections can be as large as 50% (30%) in the limiting spec- 
trum but no more than 20% (10%) for A = 1. This is not surprising since A 7^ (Qo 7^ ^qcd) ^''^duces the parton 
emission in the infrared sector and, consequently, higher-order corrections. As discussed in Sect. IIIB, the large 
and negative corrections to the color currents in the limiting spectrum lead to negative double-differential spectra, 
d'^N/dMy, at small y. Interestingly, at A 7^ 0, the infrared sensitivity is somehow weakened. As a consequence, 
d'^N/dMy is no longer negative at finite A, as illustrated in Fig. 9. Another interesting consequence is the disappear- 
ance of the infrared divergence at j/ = in the limiting spectrum, coming from the running of as', since Qq 7^ ■'^Qcn' 
as and therefore d^N/d£dy remain finite over the full momentum-space. 
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Figure 8: MLLA (left) and NMLLA (right) normalized corrections to the LO and MLLA color currents, respectively, for 
different values of A. 
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Figure 9: MLLA double differential distribution for a quark jet Ysq — 6.4 computed at A = (solid line) and A = 1 (dashed 
line) . 



C. fcj^ -distributions 



The absolute A;^ -distributions of "massive" hadrons is computed in Fig. 10 (left) for various values of A for jets 
with hardness Q = 119 GeV. As expected, as A gets larger, soft gluon emission is strongly suppressed such that 
the distribution flattens at small k± , while more hadrons are produced at large , making in turn the distributions 
harder. We also compare in Fig. 10 (right) these calculations with CDF preliminary data, all normalized to the 
log(fc^/lGeV) = —0.1 bin as before. The best description is reached in the limiting spectrum approximation, or at 
least for small values of A < 0.5. This is not too surprising since these inclusive measurements mostly involve pions. 

Predictions beyond the limiting spectrum were shown to describe very well the hump-backed shape of the inclusive 
spectra for various hadron species; in particular, the hadron-mass variation of the peak turned out to be in good 
agreement with QCD expectations (see e.g. [2]). The softening of the fc^^ -spectra with increasing hadron masses 
predicted in Fig. 10 is an observable worth to be measured, as this would provide an additional and independent 
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check of the LPHD hypothesis beyond the hmiting spectrum. This could only be achieved if the various species of 
hadrons inside a jet can be identified experimentally. Fortunately, it is likely to be the case at the LHC, where the 
ALICE [21] and CMS [22] experiments at the Large Hadron Collider have good identification capabilities at not too 
large transverse momenta. 
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Figure 10: Absolute (left) and normalized (right) inclusive fc^ -distribution beyond the limiting spectrum approximation at 
NMLLA in a jet of hardness Q = 119 GeV. 



V. 2-PARTICLE CORRELATIONS 
A. Correlators and evolution equations 

We work, like in [5], with the normalized correlators 



(7(2) n(2) 



where Gi,Qi,i = 1,2 are the inclusive spectra relative to the outgoing hadrons hi and /12, and G'-^-' , Q^^-* are the 
2-particle distributions in gluon and quark jets, respectively. The former are obtained by a single differentiation of 
the "MLLA" generating functional Z, and the latter by differentiating it twice [5] (see also the discussion introduced 
in II). Z satisfies the evolution equation described in section (2) of [5]: dZi/dlnO for the jet initiating parton A 
is expressed as an integral over z involving the DGLAP splitting functions ^^^(z) and Zb and Zc associated to the 
products B and C of the splitting process; B carries away the fraction z of the energy i? of A and C the fraction 
(1 — z) (see Fig. 11). The topology of Fig. 11 respects the exact AO constraint over the successive emission angles of 
partons (6 > 6i > 02)- In practice, suitably differentiating the master evolution equation for Za, which arises as a 
consequence of exact AO in parton cascades, yields, for the correlation functions [8] 

G(2)-GiG2 = (C3-l) Gi G2, Q(') - Qi Q2 = (C, - 1) Qi Q2, (47) 



the system of coupled evolution equations: 
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+ (Gi{z) - Qi) (q2(1 - ^) - O2) + (g2(z) - Qa) (^1(1 - z) - Qi) 

(G(2)_GiG2), = j\z^^l{z) (g(2)(z)-zG(^)) + (Gi(z)-Gi)(g2(1-z)-G2) 
■ dz^rif^liz) 2(q^^\z) - Qi{z)Q2{z)) - (g^^) - G1G2) 

■(2gi(^) -Gi) (202(1 -^)-G2) 



(48) 



+ 
+ ( 



(49) 



The derivative is taken with respect to y = Y — £ rather than with respect to InO, since it is more convenient 
when a colhnear cutoff is imposed (sec Section (2.1) of [5]). Like for the inclusive spectra, the notations have been 
hghtened to a maximum, with G^^^ standing for G^'^^z = 1) and likewise for Q^'^h The notation Xi,£i, . . . refers to 
the £i = ln(l/a;i) of the outgoing parton (hadron) i. 

B. Including NMLLA corrections 

We follow the same logic, exposed in Section II C, for the 2-particlc distributions Q^'^\ G'^^\ as the one used for the 
inclusive spectra B in Section III A. Therefore, the expansion at small a;i,a;2 is performed for ^Qi (^) (^) 
and ^Gi (^) ^G2 (^) as well as for £2^(2) xs.) and ^is^G^^) (2^, i^), similarly to Eq. (11). 



1. Quark jet 

Operating like for (16) and (17), the first (MLLA) term in the r.h.s. of (48) can be cast in the form 



Jo 



+ 



Az 



as 



$^(Z) [G(2)(z)+fQ(2)(-^_^)_Q(2) 



Cf_ 



4 7v:^°^ 



Cf 
N, 



7 CV/5_7r^ 

8 iV^ V 8 6 



Gf 



7o(GiG2)^, 



where we have plugged the DLA formula [19] 



(50) 



(51) 



18 



in the r.h.s. of (50); the terms in (51) of relative order 0(70) are neglected because their contribution provide 
corrections to (50) beyond NMLLA (see also appendix H). The second and third terms in the r.h.s. of (48) provide 
the NMLLA correction: 

dz ^ (Gr{z) - Qi) (Q2(1 -z)-Q2) = ^ dz mz) ln(l - z)j {G, - Q,)Q2,e 



Nr. 



(52) 



where the DLA expression Qi — jf-Gt + (^[lo) is used [19]; further corrections {0{%)) to this formula are here 
again dropped out because their inclusion goes beyond the present resummation logic. Likewise, we have 



/ci,^*2W(«^W -<?.)(*(! -.)-<?■) = 



Gathering (50), (52) and (53) yields 
which is written in a form similar to (18). 



Cf 

Nr. 



8 iVp I 8 ~ 6 



^2^(2) 



(53) 



(54) 



2. Gluon jet 

The structure of (49) can be worked out in the same way. The first integral term in its r.h.s. is the same as that 
in (19), such that we can simply set 



J\z ^^liz) (g^^Hz) - = £ ^7oG'('^( 
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-12^0 



2,-(2) + / 67 _ TT^ \ ^,2^(2) 



36 



The second term provides a contribution 



^GuG2t I dz^l{z)\Tiz\n{l-z) 



llTT^ 395 

h 2C(3) 

36 108 ' 



llGuGit = O(7o), 



that is beyond NMLLA and therefore dropped out here. The second line of (49) simplifies to 



dz ^nf<Pl{z) 2(g(2)(z) - Q,{z)Q2{z)) - (g^^) - G1G2) 



SNr 



7o 



X 2(g(^) - Q.Q,) - (g(^) - G1G2)] - Y8^i'o{Q''' - QiQ2)^, 
and the third one gives 

dz ^ nf^liz) (2Qi(z) - Gi) (202(1 -z)- G2) = ^7o (2Q1 - Gi) (2^2 - G2) 

-§^7o [(2Qi - Gi)02^ + (2O2 - G2)Qu . 

Gathering (55), (56), (57) and setting (see appendix H for further explanations) 



0«^G + O(7o), 



Cf f Cf 

in the subleading pieces, we obtain the NMLLA equation for the gluonic correlator 



0(^)^gG(^) + ^(^-l)G.G2+0(7o) 



(55) 



(56) 



(57) 



(58) 
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1 



3iVc 
13 n/Tfl Cf 

9 iV, A^, 



1 



1 - 

AT, 



11 + '^f^^ 

12 3iVc 



7o 



2(^(2) 



70G1G2 + 



67 
36 



_ 13 n/TijCF 
"e" ~ 18 Nr 1^ 



(59) 



The way to get the equations for the correlators Cg and Cq, to be solved iteratively, proceeds like in Section 4 and 
Appendices A and B of [5]. 



C. NMLLA correlators 

1. Gluon correlator Cg 

The differential expression for (21) reads 

Gey = ilG - ai7o' {i^e - Ml) G + a27o (V'l + ^« - /3o7o V'^) G. (60) 
Differentiating (59) with respect to £ gives the following NMLLA differential equation 

(g(2) _ GiGa)^^ = 7oG('^ - a^^l {of - /3o7oG("^) + (ai - h^)^l [(G1G2), - /3o7oGiG2] 

+ a2ll (Gf;> - /?o7oGf ^) + &27o [{GiG^)^ - Ml{G^G2)i] , (61) 

where ai, a2 are given by (24b) and (24d), and with the following coeSicients: 

Noting ^ = InG and x = InCg, the second line of (61) can be rewritten in terms of logarithmic derivatives of G and 
of Cg (see Appendix I) from which Eq. (61) is solved iteratively. Setting G^^^ = CgGiG2 in both members and making 
use of (60) leads to the analytical solution of (61), valid for arbitrary A 



C 



^ _ 1-^1-61 + il)2,t - [/3o7o]) - [QiX^ + '^2] + ^3 



9-^- : — : — 7-^. — - — z . (63) 



1 + A + (5i + 



ai (x^ + [/3o7o]) + ^2 



+ 5i 



where, like in [5] , we introduce r] = £2 — ii- Ss and ^4 are the new NMLLA corrections: 

S3{ei,£2;v) =a2fi{ii,£2;v) + b2f2ih,i2;v) = O(7o), 
S4{h,e2;v) = -a2h{iij2;v) = C>(7o), 



(64) 



and /i, /2 and /a are defined in (II) of appendix L Setting = ^4 = in (63), one recovers the exact analytical 
solution of the corresponding MLLA gluon equation (with 02 = 62 = in (61)); to derive this formula we have used 
the same method that was, for the first time, implemented in the appendix A of [5]. The other quantities and their 
order of magnitude are (see [5] ) 

X = InCg, Xi = %= O(7o'), Xy = ^ = 0{^l), (65) 

1 dG 1 dG 

= InG,, = = ^(70), ^i,v ^-G.^= ^(-^o)' (^ = (66) 



A 



1o\M2,v + i^i,y^2,i) = 0(1), (67) 



-^i = 7o X^(V'i,!/ + ^2,3/) +Xj/(V'i,€ +V'v) =C'(7o), (68) 



5. 



'2 



%\xeXy + Xiy)=0{^l). (69) 
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To evaluate (65) we consider the bare correlator: 

^ , L 1-61 + ^2,£) + [6i/?o7g] " 
^ '""['^ 1 + A + [ai/3o7o'] 

the derivatives of which are calculated numerically to eventually determine (68) and (69). 

The analytical result (63) for Cg will be numerically displayed for the limiting spectrum A = in section VD by 
using (25). For the case A ^ 0, we report the reader to [6] where it has been treated in MLLA by the steepest descent 
method. 



2. Quark correlator Cq 



The differential expression of (18) reads 



Cf 
N, 



Differentiating (59) with respect to £ gives the NMLLA differential equation 



- Po7oG| 



(70) 



(71) 



to be solved itcrativcly. Setting Q^^^ = CqQiQ2 in both members and using (70), one gets the analytical solution 
(71), valid for arbitrary A 



1,< + V'2,€ + 



M-[/3o7o]) +^3] 



Cf Gi Cf G2 ? r? 1 



A + [1 - I (Vi,^ - [/3o7o']) + 54,1] ^ + [1 - ! {i^2,e - [/3o7o1) + ^4,2] ^ § + 5i + [^2] ' 
where ^3 and 84 are the new NMLLA coefficients (02 is given by (24c)) 

4(4,4;r7) = 52./i(^i,^2;??) = 0(7o), 
Kii^uh; V) = 02/4(^1, ^2; v)- = O(7o)- 



(72) 



(73) 



Setting ^3 = 54^i = in (72), one recovers the exact analytical solution of the corresponding MLLA quark equation 
{5.2 = in (71)) that was obtained in the appendix B of [5]. We have introduced (see [5]) 



^1 = lo'^ (^e{<Pi,y + V>2,y) + (ry{(pi,e + ip2,e) =0(70), 
S2 = %'^(aeay + aey^ =0{^^), 



with ifk = InQk and a = InCg. For the numerical computation of a, we take 



In^ 1 





1- f(''J-.' + 


■•2^/ + [\/ 




Cp G, Cp G-_ 

A, Qi A, Q: 




A + 


i-Kv-M-M) 


Cf Gi , 

iVe Ql 


l-|(V'2/-[«) 


C F Gi 
Ql 



(74) 
(75) 
(76) 



(77) 



in which one uses the NMLLA expression (26) for G and Q deduced from (22) and (23), and the exact expression 
(63) {ovCg{h,y2,Ti). 

The numerical solution of (72) is given in section VD for A = 0. We make the approximation ipi k, tp^, ipy w ipy 
that is justified in Appendix E through (E2). We can therefore also use (25). The case X ^ was also dealt with at 
MLLA for a quark jet in [6] . 

Finally, taking xi = X2 in (63,72) and going to the asymptotic limit Q ^ 00 (F ^ 00), one finds the implicit overall 
normalization of these observables to be given by those of the multiplicity correlators [23] 
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Y^oo {ng{ng - 1)) 4 Y^oo {nq{nq- 1)) Nr, 
— * ; rTi — — , C„ — » ; — r7i — i 



K)2 -3' (n,)2 -'+3C^' 

for the gluon and quark jets respectively. The statement above can be easily explained; the asymptotic expressions 
of (63,72) are respectively the DLA formulae(see [8]) 

Y^oo 1 r, I \ 1 

Cg[Xx,Xi) K. 1 + T-, Cq[XY,Xi) K. 1 



and A(a;i,a;2) = 2 for x\ = X2 in the same limit. 



D. NMLLA corrections versus MLLA 



Throughout this analysis, we have consistently incorporated a set of NMLLA corrections. These were not calculated 
in the previous work [5] which was done at MLLA accuracy for A = 0. The philosophy and the basic technique are, 
however, the same (as well as in [25]). We comment below on the role of these corrections for 2-particle correlations. 
Both Ss and ^3 are dominated by their leading term, such that 

S3 « (02 + &2)(V'i/ + i^2,ef = C(7o)> '^3 ~ a2(V'i/ + V'2,£)^ = C(7o)- 

Since both 02 + 62 and 0.2 are positive and ipi increases as i Q, NMLLA corrections are expected to increase the 
MLLA solution of [5] in the limit I1+I2 ^ 0, as can be seen in (63) and (72). Thus, as found for the single-inclusive 

fcj^ -distribution, the (.xi,.T2) domain in which the two particles arc "correlated", i.e. Cg^q — 1 > 0, becomes larger 
than at MLLA. In the limit £1 +-^2 2F, the role of the new corrections is, on the contrary, expected to vanish since 
ipe^O when £^Y. 

This is indeed what appears on Figs. 12 and 14, which compare the MLLA and NMLLA solutions at the Tcvatron 
energy scale (Q = 155 GeV). While Eqs.(63) and (72) are general analytical solutions of the evolution equations at 
A 7^ 0, the numerical results displayed below are calculated at the limiting spectrum A = 0, by plugging the formula 
(25) for the inclusive spectrum into (63) and (72). The four lines in Fig. 13 show the positions in (£1,^2) space 
corresponding to the curves of Figs 12 and 14. The two upper curves of Fig. 12 correspond to line 2, its two lower 
curves to line 1; the two upper curves of Fig. 14 correspond to line 3 and its two lower curves to line 4. The 
correlations displayed in Fig. 12 and 14 appear more important in NMLLA than in MLLA. Physically, because the 
recoil of each emitting parton is better taken into account in the former approximation, less energy becomes available 
and the multiplicity of emitted particles is expected to decrease. Consequently, inside a bunch of a fewer number of 
particles, two among them get more correlated. 



E. Dependence on Aqcd 

We have tested the dependence of the gluonic correlation function Cg on A^^^ , by varying it from 150 MeV to 500 
MeV. The results are displayed on Fig. 15, as a function of £1 + £2 (left) and £1 — £2 (right). Variations are seen to 
stay below 10%. 



F. Comparison with Fong and Webber MLLA predictions 

The comparison with the predictions by Fong and Webber [24] is also instructive. Let us recall that their calculation 
is done at MLLA, yet obtained from the exact result of [5] when the two outgoing partons are taken to be close to 
the peak of the inclusive spectrum, and when the exact solution is expanded at first order in From the present 
results and that of [5], we can conclude that: 

• the convergence of the series obtained by expanding the exact MLLA result in powers of y/a^ is very bad; if one 
proceeds in this way, NMLLA corrections may be as large as MLLA, making the series meaningless; note that 
similar conclusions have been obtained in [9] when dealing with recoil effects and, more precisely, with the role 
of exact kinematics in the bounds of integrations; 
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Y=6.4, li -l2 = 1-2 Y=6.4, li -l2 = 1-2 




Figure 12: 2-particle correlations for a quark jet (left) and a gluon jet (right) as a function of £i + £2 for £1 — £2; the MLLA, 
NMLLA and Fong and Webber [24] predictions are shown as solid lines. 



• instead, in the procedure that has been adopted here, i.e. finding exact NMLLA solutions of the (approximate) 
MLLA evolution equations, NMLLA corrections turn out to be under control and their inclusion brings the 
predictions closer to Fong and Webber's. 

The present study, together with [5] , consequently stresses out the importance of dealing with exact solutions of the 
evolution equations in jet calculus. 

VI. DISCUSSION AND SUMMARY 
A. Discussion 

Energy conservation is a fundamental issue in jet calculus. While it is well known that the complete neglect of the 
recoil of the emitting parton leads to DLA (taking only into account the singular parts of the splitting functions), the 
MLLA, in which "single logarithms" are added to DLA, takes partial account of the recoil. Corrections appearing at 
higher orders in an expansion in powers of y^oj come from (i) the shifts by Inz and ln(l — z) in the arguments of the 
hadronic fragmentation functions; (ii) the non-singular terms in the splitting functions; (iii) the running of a^. Our 
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Figure 13: Positions in (^1,^2) space of Figs. 12 and 14. 



line of approach in this paper was accordingly the following: 

• we considered MLLA evolution equations as kinetic equations of QCD, and expanded their (exact) analytical 
solutions in powers of ^/al up to the order O (««)• Contributions that do not fit into such an framework are 
discarded; 

• we stuck to the logic advocated in [14, 25, 26] that, at small x and for |ln2;| ~ | ln(l — z)\ <C ln(l/a;), the 
successive corrections, MLLA, NMLLA . . . , which better and better account for energy conservation, are taken 
care of by a systematic expansion in powers of In^; and ln(l — z). 

The size of the NMLLA terms depends on the precise definition of A : a rescaling of A would change the terms 

^ '■ QCD ° QCD ° 

at this order. Systematically solving this problem would require a 2-loop calculation which has not been obtained so 
far for any multiplicity-related observable. We therefore have to consider here A^^^ as a phenomenological parameter. 
The sensitivity of our results to variations of A have been studied and found moderate (20% for inclusive k±_- 

distribution and less than 10% for correlations) when A ^ 2A 

' QCD QCD 

We left aside the question of the matching of the two definitions of the jet axis, "inclusive" direction of the energy 
flow in this work, and "exclusive" fixing from all outgoing hadrons in experiments. 

Last, hints that NMLLA corrections that has been considered here arc the dominant one can already be found 
in the work [10] where this type of NMLLA recoil effects was shown to drastically affect particle multiplicities and 
particle correlations through a factor proportional to the number of partons involved in the process. This however 
only concerns a priori 2-particle correlations. Spanning a gate between KNO phenomenon and the techniques that 
we have used here stays a challenging task which we hope to achieve in the future. 

Since calculated NMLLA corrections proved to be quite substantial, a natural question arises concerning the mag- 
nitude of higher order corrections. There, in correlation with the remarks at the end of the introduction of section 
II, it seems legitimate to consider that, since this observable is mainly sensitive to soft particles, the corrections are 
expected to be moderate. This can be different for integrated quantities like multiplicities. 



B. Summciry 

In this work, we have computed next- to- MLLA (NMLLA) corrections to the single-inclusive fc^ -distributions as well 

as 2-particlc momentum correlations inside a jet at high energy colliders. It comes as a natural extension of [4] and 
[5] in which MLLA results are provided. In particular, it exploits the same logic of using, at small energy fraction x of 
the emitted hadron, exact solutions to (approximate) evolution equations for the inclusive spectrum. The technique 
used is based on a systematic expansion in powers of -^/oj which neglects non-perturbativc effects. Nevertheless, 
it proves to be remarkably efficient to describe the preliminary measurements of (the shape of) the A;^ -differential 
inclusive cross section performed by CDF [7]. This is an indication that non-perturbative contributions play a small 
role in this observable, and concentrates in the overall normalization (LPHD hypothesis is tantamount to stating that 
in this universal factor lies the trace of the (soft and local) hadronization process). The transition from MLLA to 
next-to-MLLA enlarges considerably the domain where the computations agree with the experimental data, both in 
the transverse momentum of hadrons and in their energy fraction x. 
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Y=6.4, li + 12 = 6.4 



Y=6.4, 11 + I2 = 6.4 
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Figure 14: 2-particle correlations for a quark jet (left) and a gluon jet (right) as a function of £1 — I2; MLLA, NMLLA and 
Fong- Webber prediction. 



In our analysis, single-inclusive x-distributions as well as fc^ -spectra have been determined exactly beyond the 
limiting spectrum approximation, i. e. for arbitrary Qq 7^ A . This should in particular be relevant when dealing 
with distributions of rather massive hadrons [20]. In this respect, experimental identification of outgoing hadrons 
could provide precious additional tests of LPHD and of the physical interpretation of the infrared cutoff Qo as the 
"hadronization scale". As far as 2-particlc correlations inside a jet are concerned, future results from LHC, in addition 
to the ones of OPAL [27] and recent ones from CDF [28], are waited for to be compared with the NMLLA predictions 
presented in this study. 

The limitations of the method are in particular: 

• neglecting non-perturbative contributions may prove less justified for not so inclusive observables. In that re- 
spect, forthcoming data on 2-particle correlations from LHC promise to be very instructive. While incorporating 
some non-perturbative contributions is not excluded a priori, a systematic way to handle them is of course still 
out of reach; 

• the absolute normalization of the distributions, which involve non-perturbative effects (hadronization) is not 
predicted; 

• the calculation is performed in the small- a; limit and extrapolation to larger x may become problematic. The 
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transition to larger or from MLLA to DGLAP evolution equations, is undoubtedly also a very important 
issue. It may be tempting to proceed in this direction by going to higher orders in the expansion initiated 
in [4-6] and extended in the present work. However, the universality of MLLA evolution equations as kinetic 
equations of QCD should be cast on firmer grounds. 



C. Perspective: going to larger x 

A Taylor expansion, when used inside evolution equations, was already advocated long ago to better account for 
energy conservation [25, 26]. It appears fairly easy to realize that pushing it at higher and higher orders of Inw at 
small X inside the convolution integral (29) should play a role in it extending the domain of reliability of the solution 
to larger and larger values of x. Indeed, in (29), one integrates from u = xtou=la certain function F{\x\u — Inx). 
F is expanded at large |lna;| around |lnii| = 0, which corresponds to w = 1. If one increases x, the domain of 
integration shrinks closer and closer to its upper bound u = \. Suppose that we set a; = 1 — e. The integral becomes 
duF{\nu — ln(l — e)). Now, in the argument of F, for all u in the domain of integration \ \nu\ ~ | lnx|, such that 
a reliable expansion of F, if it exists (it depends of its radius of convergence), must involve a large number of terms. 
This is like expanding a function /(t — a) around f{—a): for \a\ ^ |lnx| ^i^lnw^lnl, a few powers of t provide a 
good approximation to f{t — a), but when a decreases, expanding f{t — a) around /(—a) uses an expansion parameter 
t of the same order of magnitude as a itself. We conclude that increasing x requires going to higher and higher orders 
in the expansion of F in powers of | Inuj. Conversely, going to higher and higher order in this expansion is expected 
to yield a solution valid in a larger and larger domain of x. 

When applied to the evolution equations themselves, and to the similar expansion in powers of (In 2) that we 
did in section II, the same kind of arguments apply, which are not unrelated with the link between MLLA and 
DGLAP evolution equations. Since NMLLA corrections to 2-particle correlations, unlike the ones for the inclusive 
fc_L distribution, are directly connected with NMLLA corrections to the evolution equations themselves, it is worth 
giving a few comments concerning this issue. 

a/ That MLLA evolution equations (4) and (5) are, at least for inclusive enough observables, valid in a much broader 
X domain than expected has been known for a long time [8]. It was furthermore noticed some years ago [15] that, 
for parton multiplicities, the exact numerical solution of MLLA evolution equations perfectly matched experimental 
results in a very large domain, and that, accordingly, the MLLA evolution equations contain more information than 
expected and the problems of finding their analytical solutions are essentially of technical nature; 

b/ at small x MLLA evolution equations are identical to DGLAP evolution equations but for a shift by Inz (z is 
the integration variable) of the variable Y ~ y + i which controls the evolution of the jet hardness [5, 8]; 

c/ for soft outgoing hadrons [x small <^ |^| = |lnx| large), this shift is negligible in the hard parton region 
(I In^l < I lnx|). However, when going to harder hadrons, that is when x grows, |^| decreases and | lnz| is no longer 
negligible. When it is so, the function to integrate is no longer safely approximated by its 0"' order expansion 
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(corresponding to Inz = 0) and higher powers of Inz are needed. This provides, in addition to the argumentation at 
the beginning of this subsection, another link between this expansion at higher orders and going to larger x; 

d/ accordingly, the Taylor expansion that we used inside MLLA evolution equations, which extends their "validity" 
to larger x, may contribute to spanning a bridge between them and DGLAP evolution equations (see for example 
[29, 30]). 

Acknowledgments: It is pleasure to thank Yu.L. Dokshitzer, I.M. Dremin and W. Ochs for illuminating discussions. 
We also thank S. Jindariani (CDF) for an invaluable exchange of information concerning CDF data. 
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Appendix A: NMLLA CORRECTIONS NEGLECTED IN THE DERIVATION OF THE APPROXIMATE 

EQUATIONS FOR THE INCLUSIVE SPECTRUM 

To get a self-contained equation for the inclusive spectrum inside a gluon jet, one needs consistently to plug in 

Q = ^[l + (ai-ai)V^]G + 0(7o'), 

Q( = ^G, + 0{^1) (Al) 

respectively, in the first and second terms of the r.h.s. of (20). Taking into account the correction proportional to ■^^ 
in (Al) would provide an extra term 

, 2nfTRCF, 



3 K. N, 



■(ai - ai)Ge, 



which adds to the r.h.s. of (21) and slightly changes the value of (24d) from 0.06 to 0.08; this number is also small, 
such that the approximation that we justify in Appendix B keeps valid. 



Appendix B: STEEPEST DESCENT EVALUATION OF (23) FOR CONSTANT 7^ 

We solve the self-contained gluon equation (23) with frozen as by performing the Mellin's transform 

The contour of integration (C) lies to the right of all poles, and ^(w, v) is the "propagator" in Mellin's space. Plugging 
(Bl) into (23) yields 



+ 7o ( ^ - «i ) y + a27o<^y 



(B2) 



The simplest way to estimate the previous Mellin's representation is by substituting the DLA saddle point loq = 70 VT 
into the MLLA (oc ai) and NMLLA (cx 02) terms. Doing so, the steepest descent evaluation of the inclusive spectrum 
at fixed as in the limit £ ^ 1 (a; <C 1) leads to 

G{1, y) « ^y^^exp (^270 - ai7o'y + «27o Vsjj) ■ (B3) 

The result, plotted on Fig. 16 together with the DLA and MLLA results (still at fixed as), shows no significant 

difference between the MLLA and NMLLA solutions. We can therefore safely use the exact MLLA solution (25) to 
compute the NMLLA inclusive fc^ -distribution. Likewise, the logarithmic derivatives = G^/G and = Gy/G 

i^i{i,y) = lo^l^^-\a2ll[jf'\ V'y (^, 2/) =70^^-0170 + ^ 0270 (B4) 

which are used to evaluate two-particle correlation, are displayed on Fig. 17 as a function of £ = Fe — y. There, again, 
the difference between MLLA and NMLLA is negligible, such that the exact MLLA expression of the single inclusive 
distribution can be taken as a good approximation in the evaluation of NMLLA two-particle correlations. 



Appendix C: SECOND DERIVATIVE OF THE SPECTRUM Gu AT A = 



The expression of the second derivative of the inclusive spectrum for a gluon jet reads 



Gu{i,y) = G{i,v){^lj + i^g,tt){i,y) 



■ + y 



+ 



/3o 



G,(£,y)- 
1 6 



sinha^s_l_i 



(CI) 



(£ + y)- 



sinh aJ^B 



Ib is the modified Bessel function of the first kind. Gu is displayed in Fig. 18 as a function of y for three values of t. 
We notice that it is negative at small values of y and gets positive at larger y. 
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Figure 16: Single inclusive spectrum at fixed Os as a function of y at Ysq = 7.5 and £ = 2.5. 
Y = 7.5 Y = 7.5 
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l = ln(1/x) l = ln(1/x) 

Figure 17: Logarithmic derivatives (left) and tpy (right) of the inclusive spectrum G{(.,Ye„) at Yq^ = 7.5. 



Appendix D: SCALING VIOLATIONS 

Varying uEQ to EQ in the argument of the DGLAP spUtting function D^^_^ in (37) yields corrections of relative 
magnitude 0{as) which are accordingly NMLLA. We need to estimate 

duuDi^{u,Eeo,uEe)= f duuDi^(u,^{u)^ (Dl) 



where 



1 / In ^ \ 



A 



Writing 

d 



\nu 
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Figure 18: Second derivative of the single inclusive spectrum as a function of y at Yeo = 7.5. 




where 



leads to 



1 1 



d\n{EQ) d^d\n{Ee) b\n{Ee)d£,' 



Inu- 



g dlTi{EQ) 
Finally, (Dl) can be approximated by 



1 - - 



1 Ini 



b ln(i;e) d^ 



I duuDi^{u,EQo,uEQ) w f duu Di^{u, EQq, EQ) 



1 1 



duu In? 



We can now estimate the order of magnitude of this correction, taking, for example, the analytic form of -D^(u) 
(non-singlet combination of quark distributions) in the w — > 1 limit [8, 31] 



with ^ = ^(li = 1) — i In (^f") ■ We need to compare 



du ii(l — u) 



with 



81 - 



4Cf 



duulnu ln(l - u) (1 - u)-i+4C^-«("=i). 



Taking, for instance, ^(^eo =6,10=3) = 0.08, which is a typical value at LEP or Tevatron energy scales, one finds 
SI /I w 0.04. When Yq Yeo- this ratio tends very fast to 0, such that the role of this correction at larger k± is 
negligible. 
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Appendix E: EXACT VERSUS APPROXIMATE NMLLA COLOR CURRENTS 



Using (26) yields the following exact (in the sense that it takes into account all subleading corrections coming from 
(26)) expression for the color currents 

(C)— *(f,y) = (C)~(£,y) + («)f(£,y)^[(ai-ai)V',,.(^,y) 
+ (ai (^ai - ai^ +0.2- {-tjjgj + ipgM) (^, u) 

+ <5(^)?(^,y)^(ai-ai)v^£(^,y) + 0(7o), (El) 



where i — g,q, and (u)* is given in [5]. The approximate expression, used in the core of the paper, only keeps 
(Cf/N,) G in (26). 

On Fig. 19, the exact and approximate color currents are shown to be in practice indistinguishable, which justifies 
the use of the latter in the core of the paper. 

Y=5.3 




L0,l = 1.5 
Approx-NMLLA, 1 = 1.5 
Exact-NMLLA, 1 = 1.5 



0.5 



Figure 19: Approximate (used in the core if the paper) and exact color currents for a gluon jet. 



We also performed the following approximation to evaluate the color current: 



Gu Ge 
Gp, G 



(E2) 



In fact, (ai — di) « 0.18 and Gu/Gi ~ Gi/G = O (70). These approximations were also made and numerically tested 
in [4, 5] (see, for example. Fig. 20 in [5]). 



Appendix F: EXPRESSION OF S{C) 



NMLLA-MLLA 



A straightforward calculation that follows from (41) gives respectively, for the gluon and quark jets, the following 
results: 



5(C) 



NMLLA-MLLA 
9 



1 e+y+x 

12.7394 -1.49751 In 

9 Ye^+X 



1 e+y+x 

356.711 -0.0369486 - - In 

9 Y0„ + X 



0.260721 - - In " ' 

9 Ye„+Xj\ Ye„ + A , 

0.377382 - ^ In : ; " ' [' ] (^g.f + Vg.w), 



9 Yeo + A 



S{C) 



NMLLA-MLLA 



1 e+y+x 

-22.6479 -0.936071 In ^ 

9 + ^ 



0.164816 -i in i±l±^) (^i±l±A^ - 
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356.711 



-0.0635496 



9 Ye„+X 

where the expression and behavior of the function (ip'^ ^ + ijjg^i. 



0.154028 - 



9 Yq^ + A 
are given in Appendix C 



Appendix G: FIXING AND VARYING ^^i^ 



Our small x calculation cannot be trusted below a certain £rau\i otherwise, as shown in Fig. 20, d?N/didy gets 
negative in the perturbative domain. We give in Table 1 values of ^min as they come out from the requirement of 



Y=4.3 




Y=4.3 




Figure 20: for a gluon jet (left) and a quark jet (right) as a function of y for three values of t 



positivity, for different values of the jet hardness (and the corresponding maximal values of y). £niin is not an intrinsic 
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Table I: Values of ^min and j/max for different values of the jet hardness 



(physical) characteristic of the system under concern (gluon or quark jet), it is only an ad- hoc parameter below which 
poor credibility can be attached to the results. One notices in Table 1 that, at a given Qj the .^min for a quark jet is 
always larger than the one for a gluon jet; this only means that our calculations can be pushed to larger x for gluons 
than for quarks without encountering problems of positivity. The question then arises whether, in calculating the 
inclusive k}_ distribution of a mixed jet, one should attach the same ^min to each of its components, which can only 
be, of course, the larger one, that is, the one of the quark component, or give to each component its proper value of 
^min as given in Table 1. The simple answer to this question comes from the fact that the two choices give, in practice, 
extremely close results. Deeper considerations on which €,„iii should be chosen are thus irrelevant. 
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For the sake of completeness, we plot in Fig. 21 the inclusive gluon fc^ -distribution at Yq^ = 6, for different values 
of both at MLLA (left) and NMLLA (right). Changing £^^-^ from 1 to 1.5 modifies the NMLLA spectrum by no 
more than 20% for log(A:^ /lGeV)=2.5. At MLLA, the dependence proves much more dramatic. Like for the variation 




Figure 21: The dependence of the inclusive gluon distribution at Ysq = 6 on i'^^^. 



with A^^^ in subsection III E, that the variations with £min seem to increase with kj_ is only an artifact due to the 
normalization at the first bin. 



Appendix H: NMLLA TERMS NEGLECTED IN VB 



The approximations we have made in (58) needs further comments; one has indeed to replace Q and Q^^^ by the 
full MLLA expressions 



Q = ^[l + (ai-ai)V£]G + 0(7o') 



(HI) 



Q1Q2 



- 1 



Cf 



1 + {hi - ai)(V'£i + i'e.2) 



1 + A 



2 + A 



O(7o) 



(H2) 



G1G2 

respectively, bi is defined in (62) and 
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(H2) was obtained in [5] and displayed later in [6]. Working out the structure of (H2) after we have inserted (HI), 
leads to 

+ ^(m - oi)(*,i + *,2)(G<^' - G1G2) + - + •('(JGiGj + 0(7j). (H3) 

As already mentioned in [6], the coefficient (61 — ai), which is color suppressed, is ~ 10~^, ~ 10~^ and ~ |. 
Thus, the whole correction is roughly ~ 10~^. This is why it is not taken into account here, which allows for analytic 
results. Introducing the terms of (H3) oc (ai — ai) in the r.h.s. of (56) provides extra terms 

• • • + - + Ml'oiG^'^ - GrG,) + ^^(ai " a^hUOiG^), 

which add to the r.h.s. of (59). They are both, in particular, color suppressed, the first one by a factor oc and 

the second one, by cx l/N^. Thus, for example, taking ■i/'^ ~ 10""'^, taking into account that 2nfTfi/3 = 1 for Uf = 3, 
the coefficient 02 defined in (24d) and which also appears in the r.h.s. of (61) would be modified to the close value 
02 ~ 0.07. Finally, since in the above 

• • • + ^Sr^^ - 1) («i - «i)7o'(GiG2). « -0.01 X 7o'(GiG2),, 



3Nc Nc \N, 

62 defined in (62) would be changed to the value 62 ~ 0.17, which only represents a 1% variation. 
The derivatives of (HI) and (H2) with respect to i are therefore respectively approximated by 

Of = f ) + f - 1) (GrG2), + Oiji), (H4) 

because the inclusion of higher order contributions (coming from the derivatives of the above O(7o) terms) in the 
non-singular parts of the equations (such as (17), (20), (52) and (57)) would yield corrections beyond the precision of 
our approach. 

Appendix I: LOGARITHMIC DERIVATIVES OF THE INCLUSIVE SPECTRUM 

The logarithmic derivatives of G, that are used in Section VC, read 

C?f = CgGiG2 ixe + i>u + M , {GiG2)e = GiG2(Vi,^ + M, 



g\} — CgGxGi 



(GiG2)« = G1G2 [(V'i,« + V'2,^)^ + + V'2,«] • 
The functions introduced in (64) and (73) are the following: 

fi{hj2;v) = + '02,£ + Xe)'^ + i^i,u + '02,« - /3o7o(V'i,£ + ^^2,^ + xd+Xu = C(7o)> 
f2{h,h;v) = {i'i,e + "4^2,1^ + ipi,u + '02,« - /3o7o(V'i,^ + i'2,e) = C(7o)> 
/3(^i,^2; rj) = 2V'i,£V2,£ + 2xiii^i,e + i^2,e) + Xu + x| - Pollxi = O(7o)> 
/4(4, ^2; v) = i^le + - 0olii^i,e = O(7o)- 
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